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EXISTENCE OF SOLUTIONS FOR DEGENERATE PARABOLIC 
EQUATIONS WITH SINGULAR TERMS 

ANDREA DALL’AGLIO, LUIGI ORSINA, AND FRANCESCO PETITTA 

Abstract. In this paper we deal with parabolic problems whose simplest model is 

' v! - A p u + B^£- = 0 in (0, T) x Q, 

u(0,x) = uo(x) in fl, 

, u(t,x) = 0 on (0,T) x <912, 

where T > 0, N > 2, p > 1, B > 0, and uo is a positive function in L°°(Q) bounded away 
from zero. 


1. Introduction 


In this paper we are concerned with homogeneous nonlinear singular initial-boundary 
value problems whose model is 

,1 Vu\p 


u — A v u + B- 


u 


= 0 , 


where 0 is a bounded open set of M jV , N > 2 , p > 1 , B <G M + and A p u = div(|Vn| p_2 Vtt) is 
the usual p-laplacian. Here and below we use the simplified notation u' in order to indicate 
the time derivative of u with respect to t. 

Singular problems of this type have been largely studied in the past also for their con¬ 
nection with the theory of non-Newtonian fluids and heat conduction in electrically active 
materials (see for instance EUQ 3 ] and references therein). 

From the mathematical point of view, the non-homogeneous elliptic case has been con¬ 
sidered in a series of papers in the last few years. Consider the equation 

, |Vu| 2 


( 1 . 1 ) 


— A u + B- 


ui 


= /; 


equipped with homogeneous Dirichlet boundary conditions. Here 7 > 0 and / is a nonneg¬ 
ative (not identically zero) function in L 1 (H). The problem is obviously singular as we ask 
the solution to vanish at the boundary of il. In [ 3 j the existence of a finite energy (i.e., in 
Rq(H)) solution to problem (11.11) has been proved if 7 < 2 and for data / locally bounded 
away from zero. 

The case of a possibly degenerate datum / has been also considered. If 7 < 1 the 
existence of a solution in was proved in [0] for general nonnegative (not identically 

zero) data, while the case 7 = 1 was faced in | 18 j provided B was small enough (we also 
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mention Ha Eg cam and the references therein for a quite complete account on the subject). 
Problems as in (11.11) with possibly changing-sign data have also been considered in [11] in 
the case 7 < 1 (see also m for further considerations concerning the strongly singular 
case). 


In the evolutive case, problems as 


( 1 . 2 ) 


p u'y J 

u( 0 , X ) = Uq{x) 
u(t, x) = 0 


in (0, T ) x n, 
in Q, 

on (0, T) x dQ, 


have been considered in the case p = 2 and 7 < 1 (see [19]). If 7 = 1 singular problems 
as (11.21) have been considered in [26, 2?j for smooth strictly positive data, while degenerate 
problems (i.e. p > 2) were studied in [28] in the one dimensional case. 

We also would like to stress that, in the degenerate case p > 2 and if B > p — 1, to find a 
solution for the model problem (11.21) is formally equivalent, through a standard Cole-Hopf 
transformation, to find a large solution to the doubly nonlinear problem 


w' = A p w m , m <-, 

p — 1 

that is a solution that blows up at the boundary of 12 (see B DU m [20] for further consid¬ 
erations on this fact). 

The aim of this paper is to study existence and nonexistence of solutions for a general 
class of singular homogeneous (i.e. / = 0) parabolic problems as (|1.2I) in the limit case 
7 = 1. We will mainly be concerned with the case p > 2. 


The paper is structured as follows: in the next section we set the main assumptions, 
we state our main result, and we introduce some preliminary tools. Section [3] is devoted 
to prove existence of a solution in the degenerate case p > 2. In Section [4] we prove 
the existence of a solution in the case p = 2 provided the size of the lower order term is 
small (e.g., B < 1 in (|1.2|) h while in Section [5] we will prove a nonexistence result in the 
complementary case B > 1. In the last section of the paper we will also give an account on 
the singular case p < 2 by showing some finite time extinction results, and by discussing 
some explicit examples of evolution. 


Notations. From now on, we will set Q = (0, T) x fl and T = (0, T) x dtl. If not otherwise 
specified, we will denote by C several constants whose value may change from line to line 
and, sometimes, on the same line. These values will only depend on the data (for instance 
C can depend on N, f 1, T, B, ck , and /3 ) but they will never depend on the indexes of the 
sequences we will often introduce. For the sake of simplicity we will often use the simplified 
notation 


/ / = / dtdx , 

>Q JQ 


when referring to integrals when no ambiguity on the variable of integration is possible. 
For fixed k > 0 we will made use of the truncation functions T\ and G & defined as 


Tfc(s) = max(— k, min(s, k )), 
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and 

G k (s) = s - T k (s) = (M - k) + sign(s). 


2. Setting of the problem and preliminary results 


We consider a parabolic differential problem of the form 
u — div a(t, x,Vu) + H(t, £, u, Vu) =0 

(2.3) < u(0,x) = uo(x) 

u(t, x) = 0 


in Q, 
in 12, 
on r, 


where 

• 12 is a bounded open set in M N , with N > 2, T is a positive number; 

• a(t, x, £) : (0,T) x!lx R^ —> R N is a Caratheodory vector-valued function such 
that 


(2.4) a(t,x,£) ■ £ > a|£| p , 

(2.5) |a(i,a:,£)| </3|£| p_1 , 


(2.6) (a(t, x, £) - a{t, x, r?)) • (£ - 77) > 0 , 

for a.e. (t,x) £ (0 ,T) x 12, for every £,77 £ R N , with £ 7 ^ r). where a, (3 are positive 
constants and p > 1 ; 

• H(t, x, s, £) : (0, T) x 12 x R + x R^ —x R is a Caratheodory function such that 

l£l p 

(2.7) 0 <H(t,x,s,0<B—, 

s 

for a.e. (t,x) £ (0 ,T) x 12, for every s > 0, £ £ R N , where B is a positive constant; 

• u o(x) is a function in L°°( 12 ) such that uq > c > 0 almost everywhere on 12 . 

In (|2.3I) . we denote by v! the partial derivative with respect to time, while Vu stands for 
the gradient with respect to the space variable x. 

Consider problem (12.31) . Here is the meaning of weak solution for such a problem. 


Definition 2.1. A weak solution to problem (12.31) is a function u in L p (0,T;Wq’ p (Q)) D 
(7(0, T; 27(12)) such that for every w CC 12 there exists c u such that u > c u > 0 in (0, T)xw; 
furthermore, we have that 



a(t, x, Vu) • Vip + 



H(t, x, u, Vu) (p 



for every 09 £ C}(\0, T) x 12), that is, for every C 1 function which vanishes in a neighborhood 
of {T} x 12 and of (0,T) x 912. 


Note that under assumption (12.71) . the function H(t,x,u,S7u) belongs to L\ oc (Q) thanks 
to the property of local positivity required on u. 


Our main result is the following: 


Theorem 2.2. If p > 2, there exists a weak solution to problem (12.31) . Moreover, if p = 2, 
there exists a solution if the constant B appearing in (|2.7|> satisfies B < a. 
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Remark 2.3. Assumption B < a is, in some sense, optimal ifp = 2: this will be the content 
of Proposition 15.11 below, in which we show how approximating problems may degenerate 
if B > a. Assumption p > 2 is also needed in this context since, as we will see in Section [G] 
finite time extinction can occur if p < 2, so that Definition 12.II should be suitably modified. 

Our strategy in order to prove Theorem 12.21 will rely on an approximation argument. 
The next subsection will introduce our approximating problems. 


2.1. The approximating problems. We consider the approximating problems 

— div a(t, x, Vw n ) + R(i, x, u n , Vu n ) = 0 in Q, 


( 2 . 8 ) 


u n (0,x) = u 0 (x) + - 


n 


1 


U n (t,X ) = - 


in 11, 
on T. 


n 


A weak solution to this problem is a function u n such that u n > j a.e. in Q, u n - 1/n € 
L p (0, T ; Wg 1,p (f2)) 0 C([0, T]; L 1 (D)) and u' n € L l (Q) + L p '(0,T; VP _1 ’ p, (D)), and such that 

(2.9) f (u n ,v) + f a(t, x, Vu n ) • Vu + ( H(t,x,u n ,Vu n )v = 0, 

Jo Jq Jq 


for any v G L p (0, T; Wo’ p (fi)) nL°°(Q). A nonnegative weak solution u n to problem (12.8p 
does exist. In fact, problem (12.81) is equivalent to 


( 2 . 10 ) 


v' - diva(f,x, Vv n ) + H(t,x,v n + -, Vu n ) = 0 in Q, 

n 

v n (0, x) = uo(x) in D, 

v n (t, x) = 0 on T, 


where v n = u n - ^. _ 

To prove that a solution of (I2.10P exists, we first extend H(t,x,s , £) to zero for s < 0, 
and, for e > 0, we consider the problem 


( 2 . 11 ) 



rj ~i / \ 

div a(t, x, Vv n e ) + £ Vn ' s) H(t,x,v n , e + -,Vv n e ) = 0 
e n 

V n ,e(0,x) = Uq (x) 


in Q, 
in Q, 


v n ,e(t,x) = 0 on T. 

A nonnegative solution v n tE to problem (12.111) exists by the results proven in pUj . Then, 
if we take a sequence of values £ | 0, one can prove that the sequence {v Uj£ } e converges 
strongly in L p { 0, T ; Wq’ p (D)) to some function v n . Then one can pass to the limit for e 0 
in the hrst two terms of (12.111) . in the sense of distributions. As far as the third term is 
concerned, we observe that, on the set {v n > 0 }, the function Is&hA converges a.e. to 1 , 
while on the set {v n = 0} (where we cannot identify the limit of w j 2ere Vu n = (1 

a.e. by Stampacchia’s result contained in [23]) the term H(t,x,v n)£ + Vu n;e ) converges 
a.e. to H(t,x,v n + Vv n ), which is zero a.e. on this set, since H(t,x , ^,0) =0 a.e. by 
assumption (12.71) . Therefore, v n is a weak solution of (12.101) . 
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2.2. Basic a priori estimates. A standard argument allows us to show that some basic 
estimates on the approximating solutions hold. We collect them in the following 

Lemma 2.4. Let p > 2, and let u n be a solution to problem (12.81) . Then, 

ll w n||LP(0,T;VULP(Q)) < C , ||«n ||l°°(Q) < C , 


and 


(2.12) f H(t, x , u n , Vu n ) < C . 

JQ 

Moreover, there exists a function u in L p (0,T;Wq’ p (LI)) such that (up to subsequences) 
u n — ^ converges to u weakly in L p (0,T]Wq ,p (LI)) and a.e. on Q. Finally, 

'S/un —> V'u, a.e. on Q. 


Proof. The proof of the first two estimates is quite standard and can be deduced for instance 
as in [5j (see also 0) using the fact that the lower order term is positive. 

In order to get (12.121) one can use ^ T e (u n — ^) as test function in (12.81) . Integrating by 
parts, dropping nonnegative terms and letting e go to zero one gets, by Fatou’s lemma 

/ H(t,x,u n , Wu n ) < 

JQ 

which implies (|2.12l) . The almost everywhere convergence of the gradients of u n is a conse¬ 
quence of (|2.12l) and of a result in [7]. 

□ 


L ( u “ w+ s) 


3. Proof of Theorem 12.21 the case p > 2 

In this section we give the proof of Theorem 12.21 in the degenerate case p > 2. We wish 

to pass to the limit in the weak formulation of (12.81) . 

A key tool in order to pass to the limit will be the following one. 

Lemma 3.1. Let p > 2, and let u n be a weak solution of problem (|2.8j) . Then, for any 
wCC fl, there exists a constant c u such that 

Un > c u > 0, in (0, T) x u, for every n in N. 

Before the proof we recall some technical tools we will use. The first one is a well known 
consequence of Gagliardo-Nirenberg inequality which is valid on any cylinder of the type 
Q = (0,T) xH with bounded Ll (see for instance [22], Lecture II). 

Lemma 3.2. Let v € L p (0,T; Wg’ p (J))) n L°°(0, T; l/(Q)), with p > 1, /3 > 1. Then 
v € L CT (Q) with a = pand 

(3.13) J Q \^ J Q \ V»\ r ■ 

Finally, we will need the following local version of a lemma by Stampacchia (see [24]). 
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Lemma 3.3. Let u:{h,r) be a function defined on [0,+oo) x [0, 1], which is nonincreasing 
in h and nondecreasing in r; suppose that there exist constants ko > 0, M, p, o > 0, and 
r) > 1 such that 

„ , Muj(k, R) v 

u(h ' T)i (h-kWl-rY ' 

for all h > k > ko and 0 < r < R < 1. Then, for every r in (0.1), there exists d > 0, given 
by 


df = 


V (p+<r) 

M2~w(fc 0 , l) 77 " 1 

(1 - r) (J 


such that 

u(d, r ) = 0 . 

In order to simplify notations, we henceforth write a(V«„) instead of a(t,x,Vu n ) and 
H(u n , Vu n ) instead of H(t, x, u n , Vu n ). 


Proof of Lemma l3.II We divide the proof into a few steps. 

Step 1. There is no loss in generality in assuming that the constant B which appears 
in ()2.7I) satisfies B > ma x(a,p — 1). 

We use v = —u~ B ^ in 112.91) . for any nonnegative ip(t,x ) £ L p (0,T;Wq’ p (Ll)) r\L°°(Q ) 
which is zero in a neighborhood of (0, T) x dfl, in order to obtain 


/ 


{u' n ,u n B tp) + B I a(Vu n ) ■ Vu n u n B L ip+ 


-B- 1 , 


IQ 


/■ 

Jq 


- I a(Vu n ) ■ ViJju n B — I H(u n , Vu n )u n B -ip = 0, 


-B„ 


'Q 


from which, taking into account the assumptions (12.41) . (12.71) and B > a, one obtains 

r T r 

./ ,.—B . 


/ 


- I I a(Vu n ) -Vilju n B < 0. 


IQ 


p—1 

Therefore, if we set u n = B ~^ p w n B+1 ~ p and 7 = , we have 


(3.14) 


p -1 

r«^)-C [ a{-w~ 1 ^Vw n )-Vil)w'^ < 0 ; 
Jo Jq 


for every nonnegative ip(t, x) £ L p ( 0, T ; ,p (fl)) D L°°(Q) which is zero in a neighborhood 
of (0, T) x dfl and for some positive constant C depending only on B and p. Observe that 
B > p — 1 and p > 2 imply that 7 > p — 1 > 1. Moreover, observe that 


B+l-p / 1 \- _i 

w n (t,x) = cn p- 1 on (0 ,T) x dLl, w n (0, x) = c yuo-\—J P =: wo n , 

for some positive constant c. In particular, since Uq is bounded away from zero, then wo n 
is bounded in L°°(f2) and the values of w n blow up on the boundary as n goes to infinity. 
We look for an a priori local bound on the L°° norm of w n . 

Step 2. Local L°° bound for w n . 
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Without loss of generality we assume that 0 € 0; we will prove that the bound holds 
true in a ball B p centered at zero of radius p with 0<r<p<R<l, then a standard 
covering argument will allow us to conclude. 

We fix k > ||u?o»i||L°°(B i {)i and we define the sets 

Ak, p (t) = {x £ B p : w n (t , x) > k} , A kyf) = {(t, x) £ (0, T) x B p : w n (t , x) > k}. 

We consider a cut-off function rj(x) £ C£°(Br) such that 

0 < rj(x ) < 1, p{x) = 1 in B r , |Vr/| < — , 

and use ip(t,x) = G k (w n (t, x))ip s (x) as test function in (13.141) . where 5 = ^2 p ^l ■ Integrat¬ 
ing between 0 and t < T, and using the assumptions (12.4p and (12.51) , we obtain 


[ (w'mK 1 Gk{w n )ip s )+ I [ \VGk(w n )\ p ip 5 
Jo Jo JA k}R (r) 


< 


C 


R — r 




0 JA k}H (r) 


We use Young’s inequality in order to absorb the term \VG k (w n )\ p 1 in the right hand side; 
using the definition of 5, we get 


< 


[ l G k {w n )ip 5 )+ [ [ \S7G k (w n )\ p ip 6 

Jo Jo Ja kr {t) 

Notice that, since 0 < tp < 1, 

|V(GfcK)/)|P < C (|VG*K)|V P + Gk{ ^fy P 


c 


(R - r)P Ja 


G k (w n ) p ip'i-p+ 1 . 


k,R 


< 


c (ivc t (w„)iv+ 


so that, observing that 5 — p = ^_ p p+1 , we hnally get 


|VGfcK)/1* < —[ 

0 JA k R (r) (R r ) P JA k 


[ (w'm'K 1 G k (lV n )p S ) + 

Jo 

Since 7 > p — 1, we can use again Young’s inequality to have 

[ {w'n,™T 1 G k (w n )<p 5 )+ f [ \\7(G k (w n )p S )\ P 

Jo Jo J A k Mr) 


G k {w n ) P <p- T-r+1 


(3.15) 


< 


c 


2(7 + 1 )T Ja k 


G k (w n y +1 ip s +- 

^ (R~ r)~/ +1 ~p 


1 

- 2(7 + 1 ) 


sup f G k (w n ) 1+1 p s dx + 
* JA kiR (t) 


C 


(7+l)p 

(R — r) 7+ 1 —p 


I Ak,R 
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Now we deal with the time derivative part. Let us define, for s > 0, 

^k(s) = f G k {a)a 1 ^ 1 da. 

Jo 

Then it is easy to check that 

**(«) > -^—G k ( s y +1 . 

7+1 

that k > ||w , o||L oo (B i {)) we obtain 

[ { w 'n^ W 7T lG k(w n )(p S ) = [ ^ k (w n (t,x))lf S > —J—- [ G k (w n ) 1+1 (f S . 
Jo Jb r 7 + 1 Jb r 


Therefore we can use (13.15[i in order to deduce 


77 [ G k {w n y +1 tp 5 < 1 sup [ G k (iv n y +1 <p s +- C (7+1)p \A ktR 

+ 1 Jb r 2(7 + 1) t J AkR 

\Ak,R\ 


7 + 1 Jb r 

and we can take the supremum over t £ (0, T) on the left in order to get 

<d k (w n )''>“ <- C 


sup [ G k (w n y +1 ip s 

t JAu R (t) 


(7+1 )p 

(R — r)~t +1 ~p 


Gathering together all these facts, and using again that ip s > , we e nd U p with the 

following estimate 

sup/ (G k {w n )^y +l + [ |V (G k (w n )tp s ) \ p < - C ( 1+ i )p \ A k,R\ ■ 

t J Ak,R(t) J Ak,R (}^ — 7 *) 7 + 1 —p 

We are now in the position to apply the Gagliardo-Nirenberg inequality (|3.13|) to the func- 
tion G k (w n )(p s , with (3 = 7 + 1; recalling that 7 = 1 on B r , we obtain 

[ c k (w n y^ <- c +1)P(N+P M w\ 1+ *■ 

J A^ r ( 7 ^ — 7 ") ( 7+1 — P )N 

Stampacchia’s procedure is now quite standard. For h > k, one obtains 


that is, 


f iV+'y+l C iV+ 7 +l iV+'y+l 

/ G k (w n ) p n > / G h {w n ) p —3— >{h- k)*>-3-\A hir \, 

J A r J Ay, r 


C\A ktR \ 1+ % 


AT +7 + 1 p( 7 +l)(iV+p) ' 

(h — k) p n (R — r) (7+i-p)iv 


Therefore, if we choose tx(h,r) = |^4/j !? .|, we can apply Lemma HOI in order to deduce that, 
for every fixed p £ (r, R), \Ah A \ = 0 if h is larger than some constant C p . It follows that 

(3.16) w n < C p , a.e. on (0, T) x B p , for every n. 

Step 3. End of the proof. Recalling that B > p — 1 and the definition of w n we use 
(13.1611 to have, a.e. on (0, T) x B p 


Un = 


p — 1 


B + 1 


B + 1 P -W n ~ B+1 ~ P > ^ ' \ y C D B+1 ~ P =C n > 0. 


p — 1 
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As we said, by means of a standard covering procedure this estimate can be proven to hold 
on any set of the form (0, T) x in, with in CC ft. □ 

3.1. Passing to the limit. In order to pass to the limit as n tends to infinity, we need 
the following result. 

Proposition 3.4. We have 

u n —» u, strongly in L p ( 0, T ; W 1,p (u))), 

for every open set in CC ft. 

Proof. The sequence {u' n } is bounded in L p, (0, T; W~ 1,p ' (ft)) + T 1 (Q). Using the Aubin- 
Simon compactness argument (see Corollary 4 in [24]) we deduce that, up to a subsequence, 

u n —> u in L P (Q ), 

for some u in 1^(0, T; W 0 1 ,p (fl)) fl L°°(Q). We will prove that, for every open set in CC ft, 

(3.17) u n -A u in L p ( 0, T; W hp (in)). 

We now introduce a classical regularization u u of the function u with respect to time 
(see TO- For every u G N, we define u v as the solution of the Cauchy problem 

1 / 

-u v + u u = u, 

< * 

„ nj/( 0 ) — uq^i/ , 

where uq )V belongs to ILq ;P (0) D L°°(fl) and satisfies 

Uo,u ~> uo strongly in L 1 (12) and *-weakly in L°°(ft) , lim — II w 0 ,^|| M /1 .p(q) = 0. 
Then one has (see US]): 

G L p (0,T-W^ P m u'„ G L p (0,T;W^ p (ft)), 

\\ u v\\l°°{Q) < IMIl°°(Q) , 

and, as v tends to infinity, 

(3.18) Uv u strongly in L p (0, T; U7 q’ p (H)). 

Let ip\(s) = se As2 (with A to be chosen later). We will denote by e(u,n) any quantity such 
that 

lim lim sup |e(u, n)| = 0 . 

^—>•+00 n—^+oo 

For 0 < (j) £ C£°(£l) we have that 

(3.19) / {u' n ,ip x (u n - u v )(j>) > e(u,n), 

Jo 

(see pH] [17]). 
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Now, using (13.191) and ip\(u n — u u )(j) as test function in (12.91) . we obtain 

/ a(Vu„) ■ V(u n - u v )y' x (u n - u v )<j> + / a(Vu n ) • V<j><p\(u n - u v ) 

JO Jo 


+ / H(u n , Wu n )(p x (u n - U v )<j> < -e{v,n). 


>Q 


Moreover, if a; CC is such that supp <p C w, since u n —* u weakly in L p (0, T; Wo’ p (fi)) 

and ip\(u n — u u ) converges to <^a(k — u v ) *-weakly in L°°(Q), we have 


'Q 


a(Vu n ) • V</> <^A(«n - U u ) = e{v, n ). 

I 

If c w is the constant given by Lemma 13.11 we have, recalling that supp^ C iv, 

I VuJ p 


'Q 


Thus, 


H(u n ,\7u n )(p\(u n - u u )(j> 


< B 


'ux(0,T) u n 


-W\{u n - U v )\<j) 


< 


B f 

— \Vu n \ p \y\{u n - U v )\(j). 
C 0J J Q 


(3.20) / a(V%) ■ V(u n - u v )ip' x (u n - u v )(j) - — [ |Vu n | p |^ A (u n - u v )\<j> < e(v, n). 

Jq c u Jq 


IQ 

Then we can write 


/ a(Vu n )-V(u n -u v )ip' x {u n -u v )(j) = / [a(Vu n ) - a(Vu u )\ ■ V(u n - u v )<p' x (u n - u^cj) 

Jq Jq 

+ / aiyuv) ■ V(u n — u u )ip' x (u n — u v )4> 

Jq 

= / [a(Va„) - a(Vu„)] • V(u n - u v )tp' x (u n - u v )<j> + e(v,n). 

Jq 

Similarly, 

/ |Vu n | P |^A(«n - u u )\(f) < a~ l / a(Vu„) • Vu n I^A^n ~ U u )\(j) 

Jq Jq 

= a -1 / [a(V«„) - a(Vu v )] • V(u n - «„) |<£>aKi - u»)\<t> 

Jq 

+ I a(yu u ) ■ V(u n -u v ) \y\{u n -u v )\<j)+ / a(Vu n ) • X7u u \ip\(u n - u v )\4> 

Jq 

1 / [o(V%) - a(Vu u )] ■ V(u n - u v ) \ip\(u n - u v )\(j) + e(i/,n). 

Jq 


'Q 


= a 


Therefore, from (13.201) we obtain 


/ [o(V%) - a(Vu„)] • V( 
Jq 


B 


u n - u u ) <p'x(u n - u v ) - \<f\{u n - u v )\ (j) < e(u,n). 

L OiC U) 
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B 1 

Choosing A large enough so that -|</?x(s)| > - for every s € M, we deduce that 

ac u 2 

/ [a(Vitn) - a(Vu u )] • V (u n - u v ) 4> < e{v, n). 

JQ 

From here it is standard (see for example ED to prove that u n — u u tends to zero strongly 
in L p (0, T; W 1,p (uj)). Recalling (|3.18|) . we thus have that (I3.17|) holds. □ 

Using Proposition 13.41 we can prove the (local) strong convergence of the lower order 
terms. 

Lemma 3.5. We have 

H(u n , Vtt n ) —> H(u, Vu), locally strongly in L 1 (Q). 

Proof. Gathering together the results of Theorem 13.41 Lemma [2.41 and Lemma 13.11 we can 
apply Lebesgue’s dominated convergence theorem to prove that 

\Vun\P |Vu|P , „ , . 

- >■ -, locally strongly m L 1 (Q). 

u n u 

It is then straightforward to conclude using (12.71) and Vitali’s theorem. □ 


Thanks to all the results proved so far we can pass to the limit in (12.81) . to have that u 
is a solution of (12.31) in the sense of Definition 12.11 thus concluding the proof of Theorem 
12.21 in the degenerate case p > 2 . 

Remark 3.6. We want to point out that since the same argument of the proof of Lemma 
13.11 applies to supersolutions of problem 12.81 Thus, in particular, one could prove the same 
result of Theorem 12.21 also for nonhomogeneous problems like 

u — div a(t, x, Vu) + H(t, x, u, Vrt) = / in Q. 

< u(0,x) = uq(x) in D, 

u(t, x) = 0 on T, 

with 0 < / in L r (0, T; L q (Q)), with 

P N i 

—I- <p, r >p , q> 1 . 

r q 


4. Proof of Theorem 12.21 the case p = 2 

As before, we start from the approximate problems (12.81) . and the key result will be the 
strong positivity of u n on (0, T) X oj, with u CC fl. 

Lemma 4.1. Assume that (12.41) . (12.51) . (12.61) and (12.71) hold with p = 2 and B < a, and let 
u n be a weak solution of problem (12.81) . Then, for any lo CC fl, there exists a constant c u 
such that 

u n > Qj > 0, in (0, T) x u, for every n in N. 
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Proof. We multiply the equation in (|2.8I) by u n e i/i, where 

e = -€( 0 , 1 ) 

a 

and 4>(t,x) G L p (0,T;W^’ p (Q)) fl L°°(Q) is a nonnegative function which is zero in a 
neighborhood of (0, T) x dfl. As in Section [3l a straightforward calculation shows that u n 
satisfies 

rT 


[ {u' n , u n e ilj)+[ ci(\7u n )-V- i/juJ > 0 . 
Jo Jq 


Then, if we define 


v n = 


i-e 


i-e 


the previous inequality can be read as 

rT 


(1 - 0) e+2 [ ( v' n ,ip)+[ a(t,x,v n , Vv n ) • VV> > 0 
Jo Jq 


Q 

where we have set 

_ e_ e 

a(t,x,s,£) = s !-» a(t, x, s 1 - 0 £). 

In other words, v n is a weak solution of the variational inequality 

(1 - 6) e+2 v' n - diva(£, x,v n , Vv n ) > 0 in Q, 

i \ 1-0 

in fl, 
on T. 


«n(0,x) = 

Mt,x) = (( 13 ^) 


1—8 


Since the vector-valued function a satisfies the standard inequalities 

a{t,x,s,£) ■ £ > a |£| 2 , 

\a(t,x,s,£)\ < P\£\ , 

(a(t, x , s, 0 - a(t, x , s, rj)) ■ - rfj > 0 , 

by the standard maximum principle (see m) we easily deduce that v n > C u > 0, in any 
u CC f 2 ; this implies the result. □ 


Once we have proved Lemma 14.11 the proof of Theorem 12.21 can be concluded as in 
Section [3l using Proposition 13.41 and Lemma 13.51 which continue to hold also for p = 2, 
with the same proof. 


5. The case p = 2 and B > a 

In this section we try to explain how nonexistence of solutions (in the sense of approx¬ 
imating sequences) may actually occur for problem (12.31) with p = 2 if the absorption 







DEGENERATE PARABOLIC EQUATIONS WITH SINGULAR TERMS 


13 


parameter B is too large. For the sake of explanation, we will restrict our attention to the 
model problem 

iVd 2 

u' -Au + B ] - -= 0 in Q, 

u 

' u(0,x) = uo(x) in fi, 

u(t, x) = 0 on T, 

for B > 1. An easy rescaling argument will then allow to deal with the case B > a for 

u 

In order to provide a simpler proof, without loss of generality, we will consider a constant 
initial datum (e.g., uo = 1); the case of a general initial datum can be obtained in a 
straightforward way by comparison (recall that uq > c > 0) and elementary changes in the 
proof. 

Our result is the following. 


Proposition 5.1. Let B > 1, and let u n be a solution of problem 

(5.21) 


u'„ — A u n + " = 0 in Q, 


V*n 

U n { 0,x) = 1 


u n (t,x ) = i 


in Q, 
on r. 


Then 


U r , 


0 uniformly on compact subsets of Q . 

Proof. We first analyze the case B = 1, readapting an idea of [20]. Let u n be a solution to 
( 15 . 211 ) and define v n = — log(u n ). Then v n solves 

v' n - Av n = 0 in Q, 

v n (0, x) = 0 in II, 

„ v n {t, x) = log(n) on T. 

On the other hand it is clear that v n = (1 — in) log(n), where w is the unique solution to 


problem 


w' = Am in Q , 
m(0, x) = 1 in ST, 


w(t, x) = 0 on T. 

By the strong maximum principle we have 

0 < w(t, x) < 1 for any (i, x) in Q, 
so that v n —> Too uniformly on compact subsets of Q. Hence, 

u n = e~ Vn —> 0 uniformly on compact subsets of Q, 


as desired. 
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If B > 1, a similar argument applies; here the right change of variable is 



which transforms problem (I5.2ip into 

v' n - Av n = 0 in Q, 


< v n (0,x) = 2=i in fi, 



The conclusion then follows as before, since v n locally uniformly on compact subsets of Q, 


and u n = [(B - l)v n \ B -i . 


□ 


6. The case p < 2. Some partial results. 


In this last section we will briefly discuss the singular case p < 2. We will show that no 
solutions in the sense of Definition 12.II can be expected as finite time extinction may occur. 
Moreover, we will also show some explicit nondegenerate examples of evolution. 

6.1. Extinction in finite time. Consider 


v! — A p u + B = 0 in Q, 




u 


( 6 . 22 ) 


ii( 0 , x) = uo(x) in D, 


u(t, x) = 0 on T 


where T > 0, N > 2. Concerning the initial datum uq we assume that uq is a nonnegative 
bounded function such that uq > c > 0 almost everywhere in Q. 

We have the following result that should be compared with finite time extinction property 
for the, formally equivalent, doubly nonlinear case (see for instance (25] and references 
therein). 

Lemma 6.1. Let p < 2 , B < p—1, uo be as above, and let u be a solution of (16.221) ; then 
there exists T* > 0 such that u(T*,x) = 0. 

Proof. Let 


w = 





w^w' — A p w = 0 in Q, 



w(t, x) = 0 


so that w solves 


w( 0, x) 


in 12 , 
on T, 
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with (5 = ■ We consider w 1 (where 7 will be chosen later) as test function in the 

previous problem, and we integrate on Q t = (0, t) x with 0 < t < T. We obtain 


w 


( t )7+/?+ 1 + / |Vu;| V 1 < C, 


' Qt 


that implies 


[ w{t) 1+fi+l + [ 
JU J Q, 


w(ty + P +1 +l \vw 2 ± ^ A \ p <c. 

'Qt 


Then using Sobolev inequality 

w(ty +0+1 + 


I 

Jn 


(7+p-l)p* \ p» 

w p I < C. 


'Qt 


Now observe that, by Holder inequality, that can be used as (7 + p — l)y >7 + /3 + 1 (for 


7 large enough), 


Hence, if we define 


[ w 

JQt 


7+/3+1 < q 


p(7+/3+l) 
(7+p-l)p* ^ P*(7+P-l) 

W P 


'Qt 


v(t) = [ 
Jc 


w 


,7+/3+l 


Qt 

we have, for some positive constants C\ and C, 

7+p—1 

rj'(t) + C\rf{t) 7 +/ 3+1 < C . 

We now remark that < 1 due to the assumptions p < 2 and B < p — 1. Therefore, 

by standard ODE analysis, there exists T* > 0 such that rj(T*) = 0; this implies the finite 
time extinction for w, hence for u. □ 

6.2. An example of nontrivial evolution. As we already said, the analysis in the case 
p < 2 is much more delicate. In the previous section we showed that in some cases finite 
time extinction may occur and so no solutions (in the sense of Dehnition 12.11) can be found 
in general. In this section we provide an explicit example of solution that shows that, in 
some cases, a solution exists at least for small times. 

Example 6.2. Let < p <2 and let B be such that 

0< s < (E—iSiE—Ll <p -i, 

p * — p 

First of all, consider the solution v to problem 


(6.23) 

where, 


-A p v = 7 : 


— t/P-L 


(7—i)(p— i)+i 


7 = 


v = 0 
p — 1 — B 

p — 1 


in fl, 
on (9fL 


> 0. 
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Problem (16.231) admits a nontrivial, positive, variational solution obtained by a standard 
mountain pass procedure as soon as 

i „ (7- 1)(P- 1) + 1 ^ * 1 

p — 1 <- <p — 1 , 

7 

and an easy calculation shows that the previous relation is equivalent to 

(p*- 2 )(p-l) 


B < 


p* — p 


Consider now the function 0 = v 1 , which solves the problem 

|V0| P 

—A„0 + -B-—— = 0 hi 

9 

0 = 0 on 312, 

and define 

u(f, x) = (1 — (2 — p)t) 53 p 0(x) . 

We can compute the problem solved by u, we have 

u' - A p u + B = (1 - (2 - p)t)^~p f-A p 0 + 

u V 0 

so that u is a solution of problem 

lVd p 

u' — A p u + -= 0 in Q, 


u 


u{ 0 , X) = 0 (x) 
u(t, x) = 0 


in 12 , 
on r, 


-0 = 0 , 


that has T* = 


2-P 


> 0 as extinction time. 
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